Dynamics of self-generated, large amplitude magnetic fields following 
high-intensity laser matter interaction 
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The dynamics of magnetic fields with amplitude of several tens of Megagauss, generated at both 
sides of a solid target irradiated with a high intensity (~ lO^^W/cm^) picosecond laser pulse, has 
been spatially and temporally resolved using a proton imaging technique. The amplitude of the 
magnetic fields is sufficiently large to have a constraining effect on the radial expansion of the 
plasma sheath at the target surfaces. These results, supported by numerical simulations and simple 
analytical modeling, may have implications for ion acceleration driven by the plasma sheath at the 
rear side of the target as well as for the laboratory study of self-coUimated high-energy plasma jets. 

PACS numbers: 52.25.Xz, 52.38.Fz, 52.70.Nc 



The generation of magnetic fields in plasmas is a phe- 
nomenon of great relevance for a wide range of physi- 
cal scenarios and it has been studied in laser-produced 
plasmas since the introduction of high-power lasers, with 
particular emphasis to their role in Inertial Confinement 
Fusion In this scenario, involving nanosecond laser 
pulses of intensity 1^ ~ lO^** — 10^^ W/cm^, magnetic 
field generation is generally described by hydrodynamic 
modeling and it has been recently characterized by 
temporally and spatially resolved measurements In 
the case of shorter (< ps) and more intense {1^ > 10^^ 
W/cm^) laser pulses, the production of large currents 
of high energy "hot" electrons requires a more complex 
modeling Q and very intense magnetic field may be gen- 
erated by, e.g., Weibel-like instabilities ^ or electron re- 
circulation at the plasma boundary (fountain effect) 0. 
In this latter case, magnetic fields may influence the emis- 
sion of multi-MeV proton 0, Q and positron [l^ beams. 
Magnetic fields generated in relativistic laser-produced 
plasmas are also of central importance in reproducing 
conditions resembling large-scale astrophysical processes 
[m on a laboratory scale such as the self-coUimation of 
relativistic leptonic jets [l3 or the upstream-downstream 
mixing in supernova remnant shocks (Tsj . 

Previous experimental work has detected effects in- 
duced by such magnetic fields on external optical beams 
14 1 or on the polarization of self-generated harmonics 
15|. However, these measurements suffered of limita- 
tions, in terms either of the range of plasma density ac- 
cessible or of spatial and temporal resolution (isj . 
Furthermore, only the fields generated at the front (laser- 
irradiated) side of the target have been investigated; the 
magnetic field generation at the rear side, where proton 



acceleration in the expanding fast electron sheath takes 
place [l6j . is thus yet to be experimentally characterized. 

In this Letter we report on simultaneous measurements 
of the magnetic fields generated at the front and rear side 
of a solid target irradiated by a short and intense laser 
pulse, using a spatially and temporally resolved proton 
imaging technique [l7[. Clear evidence is given of the 
generation of toroidal magnetic fields (maximum ampli- 
tude of ~ 50 MegaGauss) that decay in time on a pi- 
cosecond time scale. Their spatial distribution and am- 
plitude is consistent with the recirculation of the laser- 
accelerated electrons around the target and they are suf- 
ficiently intense to confine the radial plasma expansion. 

The experiment was carried out at the Vulcan laser 
system in the Rutherford Appleton Laboratory [l^ us- 
ing two laser beams both with a central wavelength 
Xl = 1.05/^m, energy El — 50 J and duration tl = Ips. 
Both beams were preceded by a lower intensity plateau 
{Ip ?a lO^^W/cm^, duration of w 300 ps), due to ampli- 
fied spontaneous emission. The first laser beam (CPAi) 
was focussed, down to a focal spot with radius ~ 
10/im, to a peak intensity of II ~ lO^^W/cm^ (dimen- 
sionless intensity oq ~ 2.7) onto a d = lO^m thick alu- 
minum foil, with a 45° angle of incidence. Hydrodynamic 
simulations [l9j indicate that the laser pre-pulse ablates 
a sub-micron layer of the aluminum foil that expands, in 
a plasma state, with an exponentially decreasing density 
profile (scale-length of 2.5 microns) and an electron tem- 
perature of «100 eV. The second laser beam (CPA2) was 
focussed onto a 20yum thick gold foil to generate, via Tar- 
get Normal Sheath Acceleration (TNSA) 0], a proton 
beam with a Boltzmann-like spectrum (temperature of 
Tp = 3.0 ± 0.2 MeV and cut-off energy of Ep « 20MeV, 




Figure 1. a.-f. Proton imaging snapshots in the direct configuration (sketched in the frame h.) and typical image of the 
unperturbed proton beam (g.)- i--k. Additional set of images, still obtained in the direct configuration, from a different 
laser shot. l.-n. Proton imaging snapshots in the reverse configuration (sketched in the frame p.) and typical image of the 
unperturbed proton beam (o.). In all images, the spatial scale refers to the interaction plane and time is relative to the arrival 
of the peak of the CPAi pulse on target. 



see Figs. [T]g and[TJo for its spatial distribution). This 
beam was used as a charged-particle probe and was 
recorded onto a stack of calibrated RadioChromic Films 
(RCFs) [l^l giving a point-like projection of the interac- 
tion with a geometrical magnification M ^ (l + L)/ 1 fa 11 
with I ?a 3 mm and L w 3cm the distances between the 
Au and Al foils and between the Al foil and the RCF 
respectively, as sketched in Figs. [T]h. and[T]p. 

Previous theoretical modeling of high-intensity irradi- 
ation of thin solid targets @ indicates the generation 
of magnetic fields having a toroidal structure with az- 
imuthal symmetry and field lines parallel to the target 
surface. The intense electric fields generated in the ex- 
panding plasma sheath are almost normal to the target 
surface and have been directly detected using a probe 
proton beam parallel to the surface [2l|. In the present 
experiment, in order to maximize the probe protons de- 
flections due to magnetic fields over those due to elec- 
tric fields, the propagation axis of the probe beam was 
normal to the target surface. In order to ascertain the 
magnetic nature of the deflecting fields, two configura- 
tions were adopted: the proton beam first encountered 
either the rear, un- irradiated side of the target {direct 
configuration. Fig. [TJh) or the front, irradiated side {re- 
verse configuration. Fig. [Hp). For a given polarity of the 
magnetic field distribution, the two configurations should 
induce opposite deflections. 

A typical set of RCF images, obtained in the direct 
configuration, are displayed in Figs. [T] (scries a.-f. and i.- 
k.). All images depict the presence of two main features: 
an outer ring (radius of the order of 100 - 200 /im) and 
an inner dot (radius of the order of 30-40 /im) of proton 
accumulation. The outer ring is seen to slowly expand 
in time whilst roughly preserving its amplitude whereas 
the inner dot becomes weaker and eventually disappears 
as time progresses. The reverse configuration (Figs. [1] 
l.-n.) induces an inverse deflection pattern, the central 



dot being much darker than the outer ring. In the di- 
rect configuration, the rear magnetic field focusses the 
protons propagating near the axis whereas the front field 
enhances the divergence of the protons propagating at a 
wider angle (see Fig. [5]a). In the radiographs, this leads 
to the formation of an accumulation dot produced by the 
rear field and of an outer ring produced by the front field. 
In the reverse configuration, the same structures should 
be expected, yet with a different relative amplitude (Fig. 
[2]b), in qualitative agreement with Fig. [TJ 

For a quantitative analysis of the magnetic field. Par- 
ticle Tracing (PT) [22| simulations have been performed. 
The code follows the propagation of a proton beam, as 
employed in the experiment, through a prescribed mag- 
netic field distribution and it includes the response of 
an RCF detector providing synthetic RCF images. The 
magnetic field is assumed, at both the front and the 
rear side of the target, as having a cylindrically symmet- 
ric toroidal distribution localized at the target surfaces, 
analogous to the model reported in Ref. In this 

model, the magnetic field at each of the two surfaces is 
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Figure 2. Sketch of the proton deflections induced by the 
magnetic fields in the direct (a.) and reverse (b.) configura- 
tion. 
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described via the following parameters: B (field ampli- 
tude), R (radius of the field distribution), and L^. and 
Lz (field scalelength in directions parallel and perpen- 
dicular to the target surface, respectively). These pa- 
rameters were varied independently until a satisfactory 
match was reached between experimental and synthetic 
data, in terms of position and amplitude of the optical 
density peaks observed on RCF in correspondence to the 
ring and dot proton accumulation regions (see Fig. [3]). 
PT simulations indicate that variations in B and R at 
both surfaces affect strongly peak amplitude and radius, 
respectively, while the field scale-lengths have only a sec- 
ondary effect on the synthetic RCF profile. By keeping 
Lz = Lr = 10 /im (consistent with PIC simulations, see 
later) and varying independently B and R at front and 
rear, profile matching as shown in Fig. [3] could be ob- 
tained, in both the direct and reverse configuration. The 
matching was considered satisfactory whenever the po- 
sition and peak value of the simulated optical density 
maxima reproduced those of the experiment within ± 5 
/im and 5%, respectively. These values are of the order of 
the intrinsic spatial resolution of the proton backlighting 
p7| and of a typical small-scale non-homogeneity of the 
proton density on a beam cross section. 

By iteratively applying this method to each RCF layer, 
it has been possible to simultaneously extract the tempo- 
ral evolution of the amplitude and radius of the magnetic 
field at each side of the target (see Fig. 0]). In correspon- 
dence to the falling edge of the laser pulse, both fields 
are seen to increase their radius and decrease their am- 
plitude in time. After the laser irradiation (t > 1 ps in 
Fig. 2]) the front field is seen to rapidly drop down to an 
almost constant amplitude of Bp « 10 - 12 MG while its 
radius increases up to an approximately constant value 
of 90 /im. On the other hand, the rear field exponen- 
tially decreases in amplitude with a typical time-scale of 




Figure 3. Comparison between the experimental proton dose 
deposition as recorded by the RCF and that resulting from 
PT simulations assuming magnetic field distributions as the 
ones depicted in Fig. [5] Spatial scales refer to the interaction 
plane. Frames a. and b. refer to the direct configuration 
whereas frames c. and d. to the reverse one. 
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Figure 4. Temporal evolution of the magnetic field amplitude 
and width at the rear (a.) and front (b.) side of the target as 
extracted from matching PT simulations, compared with the 
PIC results. Time refers to the arrival of the peak of CPAi 
on target. 



the order of 7 ps. Meanwhile, its radius increases in time 
with decreasing radial velocity. The longer persistence of 
magnetic fields at the front surface might be related to 
the presence of the underdense prc-plasma which is able 
to better support the magnetic field lines [23j . 

The presence of McgaGauss magnetic fields of opposite 
polarity at both sides of the target, is also supported by 
two-dimensional (2D) Particle-In-Cell (PIC) simulations. 
A density profile composed by an exponential ramp re- 
producing the above mentioned preplasma, followed by a 
plasma bulk with electron density = 40nc and charge- 
to-mass ratio Z/A — 9/26 is assumed. The laser pulse 
has a Gaussian transverse profile (FWHM = 5 /im) , a du- 
ration of 250 T (with the laser period T = 3.3 fs for 1 /im 
wavelength), dimensionless peak amplitude oq = 2.7 on 
axis and it is incident at 20° with respect to the nor- 
mal of the target surface. The simulation ran up to 
600r ~ 2 ps, in order to overlap with the proton ra- 
diographs at earlier times. FiglS] shows the generation 
of magnetic fields at the front and rear surfaces of the 
target, both having an approximately antisymmetrical 
distribution with respect to the axis, and polarity oppo- 
site to each other. Near the peak of the laser pulse the 
fields reach a maximum amplitude of the order of 50 MG 
(Fig. 2. a) and, during the rise of the laser pulse, they 
propagate in the transverse (y) direction with a constant 
velocity of ~ 2.7 x 10^ m/s (consistently with the sce- 
nario experimentally investigated in Ref . j23| ) while, as 
the laser intensity falls down, the field distribution drifts 
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at a much lower velocity of 2 x 10^ m/s until stagnation 
is reached (t > 10 ps). Within the intrinsic approxima- 
tions that a PIC model unavoidably introduces (such as 
two-dimensional geometry and non-collisionality) a fair 
agreement is found between the experimental and nu- 
merical results (see Fig. |4]). The PIC simulations also 
indicate that the electrostatic fields at the target's sur- 
faces (not shown for brevity) are almost normal to the 
original surface (i.e. parallel to the main axis of propa- 
gation of the probing proton beam), in agreement with 
reported experimental observations |2l| , and allow to es- 
timate their intensity. The inclusion of such electrostatic 
fields does not affect significantly the PT images and can 
thus be neglected. 

In principle, strong magnetic fields may also be gen- 
erated inside the target due to resistive return currents 
which must balance the fast electron flow. In order to 
evaluate these fields, 3D simulations of the propagation 
of an electron beam through aluminium at an initial tem- 
perature of 1 eV were performed, using the code ZEPHY- 
ROS [1^. The simulation assumed suitable parameters 
for the electron beam (electron energy of 0.6 MeV, beam 
density of 4 x 10^°cm~'^, initial radius rs ~ 8.5 /im. di- 
vergence 9d — 25° [2^) and a background-temperature 
dependent target resistivity [l^l . Simulation results show 
the growth of small-scale filaments [2^ with magnetic 
fields amplitude up to ~ 40 MG and a characteristic spa- 
tial scale of 4 fj,m (Fig. 2.e) which is below the proton 
imaging resolution |17| . Indeed, including magnetic fields 
in the target bulk with these simulated amplitude and 
spatial distribution induces in PT images the superposi- 
tion of random fluctuations with amplitude below 5%. In 
these specific experimental conditions, the proton deflec- 
tions are thus predominantly induced by the magnetic 
fields generated at the surfaces of the target. 

It is of particular interest to analyze in more detail 
the field dynamics at the rear side of the target, which 
are much less explored, and are of direct relevance to 
the sheath ion acceleration. Here, the intense magnetic 
fields are generated by hot electron currents which, when 
crossing the rear surface, can not be balanced anymore 
by a countcrpropagating return current. The expected 
temperature of hot electrons produced in the interaction 
is Th — 0.6 MeV, corresponding to a velocity Vh — 0.9c 
and to a relativistic factor 7 ~ 2; the hot electron den- 
sity may be thus roughly estimated by a balance of en- 
ergy fluxes, JIl = rihVhTh yielding, for an absorption 
fraction / « 0.1 [11], n,, ~ 4 x lO^o cm^^. The total 
current due to hot electrons flowing through the tar- 
get may be estimated as Ih = erihVhS ~ 5 x 10^ A 
where S = 7rr| = 3 x 10~® cm~^ is the area of the 
laser focal spot. The generation of a magnetic field re- 
quires a diverging electron flow since, if it were colli- 
mated, the large back-holding electric field E resulting 
from charge displacement would cause an equal and anti- 
parallel displacement current J^; = eo(?tE of rcfluxing 
electrons. In this case the source term for the mag- 
netic field would exactly vanish. A divergent flow al- 



lows part of the current to flow in the radial direction 
forming loops which fall back to the target where a sur- 
face return current may close the circuit. We developed 
[2^ a simple geometrical model of such "fountain effect" 
to estimate the peak magnetic field as B^^x — ctOdBo 
where Bq = /io-^/i/(27rrQ), 9d — 25° is the divergence of 
the flow. To ~ 15 /im the radius of the electron emit- 
ting area and a ~ 8Th/{eEro) with E the typical value 
of the electric field. By estimating E ~ lO^'^ V m^^, 
as suggested both by the PIC simulations and by pro- 
ton emission data in similar conditions [2l| . we obtain 
a peak value of B„ax — 70 MG, in fair agreement with 
the experimental results. The corresponding value of the 
Larmor radius mejVh/eB ~ 3-Bj^^ /zm (where Bio is the 
field in units of 10 MG) is small enough to indicate that 
electrons would be strongly magnetized in the regions of 
peak field. In these conditions, the magnetic field lines 
are frozen in the electron fluid, and the magnetic forces 
tend to confine the plasma. This tendency may account 
for the "negative acceleration" in the torus radius ver- 
sus time cPvp^^/dt'^ < observed in Figd) The magnetic 
energy density (u™ = B'^/2^o = 4 x lO^^Bf^ J cm~^) be- 
comes in fact comparable to the maximum plasma ther- 
mal energy density expected at the peak of the laser pulse 




Figure 5. a.-c: PIC simulation results. The frames show 
the transverse magnetic field (Bx) a.t t — 442T in units of 
-Bo = rUeLJc/e — 107 MG (a.), and a zoom of the front (b.) 
and rear (c.) regions corresponding to the dashed rectangles 
in a. ai t = 442, 492 and 592T. The black arrows indicate 
the position of the field maxima, d.-e.: transport simulation 
results. The frames show the hot electron density inside the 
target in logarithmic scale and units of cm~'^ (d.) and the 
related magnetic field distribution in units of MG (e.), both 
at t = 1 ps after the peak of the pulse. 
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(ut = rihTh ~ 4 X 10^ J cm~^). Since the typical times 
for coUisional dissipation and magnetic diffusion (of the 
order of tens of ns) are much longer than the time scales 
of the observation, the magnetic field value should mainly 
decay in time because of the sheath expansion and thus 
be roughly inversely proportional to the square of the 
observed radius, in agreement with the observations. 

In conclusion, temporally and spatially resolved pro- 
ton imaging indicates the generation of toroidal magnetic 
fields having tens of MegaGauss strength on both sides of 
a foil irradiated by an intense laser pulse. The magnetic 
fields are strong enough to effectively confine the radial 
expansion of the plasma region where they are generated 
thus possibly affecting the ion acceleration in the expand- 
ing sheath. Morevoer, the self-confining effect qualita- 



tively resembles the coUimation of leptonic astrophysical 
jets, suggesting that the present framework is suitable to 
investigate similar mechanisms in down-scaled laboratory 
experiments. 
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